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GASPS — A Herschel Survey of Gas and Dust in Protoplanetary Disks: 

Summary and Initial Statistics 
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I. Tilling, 26 J. M. Torrelles, 27 B. Vandenbusche, 28 S. Vicente, 2 
G. J. White, 29 - 30 and P. Woitke 31 
Received 2012 December 31; accepted 2013 April 01; published 2013 May 31 

ABSTRACT. We describe a large-scale far-infrared line and continuum survey of protoplanetary disk through to 
young debris disk systems carried out using the ACS instrument on the Herschel Space Observatory. This Open 
Time Key program, known as GASPS (Gas Survey of Protoplanetary Systems), targeted ~250 young stars in narrow 
wavelength regions covering the [01] fine structure line at 63 /im the brightest far-infrared line in such objects. A 
subset of the brightest targets were also surveyed in [01] 145 /im, [CII] at 157 /im, as well as several transitions of 
H 2 0 and high-excitation CO lines at selected wavelengths between 78 and 180 /im. Additionally, GASPS included 
continuum photometry at 70, 100 and 160 /n n, around the peak of the dust emission. The targets were SED Class II- 
III T Tauri stars and debris disks from seven nearby young associations, along with a comparable sample of isolated 
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Herbig AeBe stars. The aim was to study the global gas and dust content in a wide sample of circumstellar disks, 
combining the results with models in a systematic way. In this overview paper we review the scientific aims, target 
selection and observing strategy of the program. We summarise some of the initial results, showing line identifica- 
tions, listing the detections, and giving a first statistical study of line detectability. The [01] line at 63 /j,m was the 
brightest line seen in almost all objects, by a factor of ~10. Overall [OI]63 /im detection rates were 49%, with 100% 
of HAeBe stars and 43% of T Tauri stars detected. A comparison with published disk dust masses (derived mainly 
from sub-mm continuum, assuming standard values of the mm mass opacity) shows a dust mass threshold for [01] 
63 /im detection of ~10 -5 M 0 . Normalising to a distance of 140 pc, 84% of objects with dust masses >10 -5 M 0 
can be detected in this line in the present survey; 32% of those of mass 10 _( -10 -5 M 0 , and only a very small 
number of unusual objects with lower masses can be detected. This is consistent with models with a moderate UV 
excess and disk flaring. For a given disk mass, [OI] detectability is lower for M stars compared with earlier spectral 
types. Both the continuum and line emission was, in most systems, spatially and spectrally unresolved and centred 
on the star, suggesting that emission in most cases was from the disk. Approximately 10 objects showed resolved 
emission, most likely from outflows. In the GASPS sample, [OI] detection rates in T Tauri associations in the 
0.3-4 Myr age range were ~50%. For each association in the 5-20 Myr age range, ~2 stars remain detectable 
in [OI]63 /tm, and no systems were detected in associations with age >20 Myr. Comparing with the total number 
of young stars in each association, and assuming a ISM-like gas/dust ratio, this indicates that ~18% of stars retain a 
gas-rich disk of total mass ~1 M Jupiter for 1-4 Myr, 1-7% keep such disks for 5-10 Myr, but none are detected 
beyond 10-20 Myr. The brightest [OI] objects from GASPS were also observed in [OI]145 /im, [CII] 157 /an and 
CO J = 18 — 17, with detection rates of 2(M10%. Detection of the [CII] line was not correlated with disk mass, 
suggesting it arises more commonly from a compact remnant envelope. 


1. INTRODUCTION 

One of the most significant astronomical discoveries of the 
past decade has been the realisation that roughly 20% of main- 
sequence FGK stars harbor planets (Fischer & Valenti 2005; 
Borucki et al. 2011). Moreover, at least 16% of FGK main- 
sequence stars are found to have a debris disk more massive 
than the dust in our own Solar System, indicating an unseen 
population of colliding planetesimals (Trilling et al. 2008; Wyatt 
2008). These two independent results imply that the planet and 
planetesimal formation process is common and robust, and can 
lead to a wide diversity of systems. However, it is not clear how 
young gas-rich disks — where planet formation is either still oc- 
curring or has recently completed — evolve into mature plane- 
tary and/or debris disk systems. 

Both debris disks and main sequence planetary systems are 
gas-poor. Debris systems, composed of grains in a collisional 
cascade, have dust masses of <10 - ' M 0 , although the mass of 
planetesimals — thought to be the starting point of the cascade — 
may be 10 -4 M 0 or more (Wyatt & Dent 2002). In most cases, 
no molecular gas is detected (Dent et al. 1995); however, in a 
very few nearby young debris systems such as /3 Pic, a small 
mass of mostly atomic gas is seen (Lagrange et al. 1998; 
Olofsson et al. 2001; Roberge et al. 2006). Possible gas forma- 
tion mechanisms in such systems are secondary release during 
grain-grain collisions (Thebault & Augereau 2005), photode- 
sorption from dust grains (Chen et al. 2007), or sublimation 
from comets (Zuckerman & Song 2012). 

By contrast, the material around young, pre-main-sequence 
(but optically-visible) stars is gas-rich. Such so-called proto- 
planetary disks are found towards ~10% of stars aged 5 Myr, 


and at least 80% of stars aged <1 Myr (Haisch et al. 2001). 
Similar in size to debris disks, then - dust masses are typically 
10 _ l -10 _i Mq. The assumption normally made is that 99% 
of the disk mass is gas, the same as that of the natant interstellar 
cloud, leading to these disks being described as ‘primordial’. 
Their total masses would then be similar to that of the minimum 
mass Solar Nebula (Williams & Cieza 2011). In the even youn- 
ger, so-called ‘protostellar’ stage (typically <0.1 Myr), disk 
masses may be still larger, approaching that of the protostar it- 
self. Systems at this phase are usually optically obscured, as the 
remnant cloud and infalling envelope have not yet dispersed. 

Dissipation of the primordial disk gas limits the timescale for 
giant planet formation, affects the dynamics of planetary bodies 
of all sizes during their formation, and determines the final 
architecture and constitution of planetary systems. No planet 
formation will take place without gas to damp the particle ve- 
locities. The methods for removal of the gas and dust compo- 
nents are generally different. Photoevaporation is thought to be 
important in gas dissipation (Gorti & Hollenbach 2004, 2009), 
and bipolar outflow jets may also play a role. Molecular species 
will also be depleted in regions exposed to the photodestructive 
effects of UV (Kamp & Sammar 2004). Selective removal of the 
dust can be caused by the interaction with a planet (Rice et al. 
2006), ice or refractory grain sublimation (Thi et al. 2005), or by 
grain growth and settling (Dullemond & Dominik 2005). As 
noted by these authors, these effects can be very rapid, occurring 
on timescales ~10 4 yr — significantly shorter than the disk ages. 

Statistical studies show that the presence of lupiter-mass pla- 
nets in mature systems is strongly influenced by stellar metal- 
licity, mass, and binary companions (Fischer & Valenti 2005; 
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Johnson et al. 2007). Do the stellar parameters also affect disks? 
Age clearly affects the fraction of stars with primordial disks 
(Haisch et al. 2001) and, on a longer timescale, debris disks 
(Carpenter et al. 2009). Disks may be affected by binary com- 
panions in debris systems (Trilling et al. 2007) but not substan- 
tially by stellar metallicity in debris or protoplanetary systems 
(Greaves et al. 2006; D’Orazi et al. 2011; Maldonado et al. 
2012). And although there appears to be no direct correlation 
between debris disks and planets (Bryden et al. 2009), there 
are clearly some systems which have both, and where the planet 
creates a gap or affects the disk shape (Kalas et al. 2008; 
Thalmann et al. 2010; Hughes et al. 2011). At present, there 
is no clear observational evidence that primordial disk lifetimes 
are significantly affected by the stellar mass (Boissier et al. 
2011; Ercolano et al. 2011). 

To study disk evolution and look for general trends, many 
large and unbiased dust continuum surveys have been carried 
out, e.g., in the near-infrared (Kenyon & Hartmann 1995; 
Haisch et al. 2001), the mid- to far-infrared (most complete 
out to a wavelength of 70 /iin in the Spitzer projects FEPS 
[Hillenbrand et al. 2008] and c2d [Evans et al. 2009]), and 
in the sub-mm (Andrews & Williams 2005). They indicate pri- 
mordial disk lifetimes of a few Myr, but there is a broad distri- 
bution of dust mass at any particular age, with notable outliers. 
For example, although the fraction of disk-bearing stars in the 
~1 Myr-old Taurus star-forming region is as high as 75%, a sig- 
nificant minority of its stars have no detectable dust excess 
(Luhman et al. 2010). 

Dust emission is ostensibly easy to inteipret, as the normally 
optically-thin sub-mm continuum can be used to directly esti- 
mate the dust mass, M d , by employing a mass opacity, and 
emissivity power law, /?. Typical values adopted in the literature 
are k v = 1.7 gem -2 and (3 = 1.0 (Andrews & Williams 2005). 
However, k„ depends on the grain size distribution (D’Alessio 
et al. 2006), and most of the solid body mass may be in large 
grains contributing little to the observed flux (Wyatt & Dent 
2002; Krivov et al. 2008). One option is to define a dust mass 
which only includes solid material smaller than 1 mm (Thi et al. 
2010); this is reasonably consistent with the standard literature 
value of Deriving the total disk mass requires an assumption 
of the gas-to-dust ratio; normally the interstellar medium value 
of 100 is used, but it is unclear whether this value is maintained 
in disks (and it is certainly not valid for debris-dominated sys- 
tems). Throughout most of this current work, we have quoted 
disk masses in terms of the dust mass, M d , allowing compar- 
isons to be made independently of the gas/dust ratio. 

Although gas dominates the mass (at least for protoplanetary 
disks), emission lines are generally more difficult than the con- 
tinuum to both observe and interpret. The bulk of the gas is in 
H 2 which has no dipolar moment. The observed intensity is 
affected by abundance variations due to complex chemistry, mo- 
lecular photodissociation, or freezeout in the cool disk mid- 
plane, as well as high optical depths and uncertain excitation 


processes. Unlike the continuum, a more limited number of 
gas surveys have been carried out. Mid-IR studies of H 2 0 show 
emission from the inner 10 AU in many T Tauri but few HAeBe 
stars (Pontoppidan et al. 2010). Spitzer surveys of [Nell] at 
12.8 /im show warm gas in many systems, thought to arise from 
winds from the disk surface at radii up to ~10 AU (Lahuis et al. 
2007; Pascucci & Sterzik 2009). At sub-mm wavelengths, lim- 
ited surveys of low-J rotational lines of CO have been carried 
out, both with single-dish telescopes (Zuckerman et al. 1995; 
Dent et al. 2005), and with interferometers (Oberg et al. 2010). 
Most of this emission arises from 30 to 300 AU radii. In the 
far-IR, scans of a few bright embedded Class 0 YSOs (e.g., 
L1448-mm) and massive young stars (e.g., Orion-KL) with 
the long-wavelength spectrometer (LWS) on ISO showed rich 
spectra, including fine structure lines of [OI] and [CII], and 
many transitions of H 2 0 and CO (Benedettini et al. 2002; 
van Dishoeck 2004). Class I YSOs also show similar lines, al- 
beit fainter than the Class 0s. Among less embedded systems, 
the bright ‘prototypical’ Class I— II object T Tauri 32 also has 
many FIR lines (Spinoglio et al. 2000). However, more typical 
optically-visible Class II — III objects were not detected in FIR 
lines due to the relatively low sensitivity and large beams. 

Both gas and dust observations suffer from problems in inter- 
pretation, and both are needed for the best understanding of disks. 
The motivation for GASPS was to conduct a relatively large, sys- 
tematic study of gas and dust in the far-infrared, utilising the 
sensitivity improvements available with the Herschel Space Ob- 
servatory. 33 The survey covers a broad sample of optically-visible 
young systems, from Class II gas-rich protoplanetary disks 
through to Class III objects and gas-poor debris disks. It focusses 
on the brightest lines and the FIR peak of the continuum emis- 
sion, and is complemented by data at other wavelengths. In this 
article, we describe the survey and observing techniques (§ 2), 
discuss the target selection criteria, and give the complete target 
list (§ 3 and the Appendix A). Section 4 summarises the origins of 
FIR line emission in these objects, and outlines the modeling 
used by GASPS. In § 5 we give an overview of the results, with 
the lines detected and the detection statistics. 

2. THE GASPS SURVEY 

GASPS (Gas Survey of Protoplanetary Systems) uses the 
photodetector Array Camera & Spectrophotometer (PACS) 
(Poglitsch et al. 2010) on the Herschel spacecraft (Pilbratt et al. 
2010) to study a predefined set of the brightest lines and dust 
continuum in the far-infrared from a relatively large sample of 
targets. The aim was to allow the detection of gas in systems 


“There is some discussion as to whether T Tauri should actually be classified 
as a Class I YSO with a massive envelope: typically, the prototypical object in a 
class actually turns out to be rather unusual! 

33 Herschel is an ESA space observatory with science instruments provided by 
European-led Principal Investigator consortia and with important articipation 
from NASA. 
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with a disk mass limit similar to, and possibly lower than, existing 
sub-mm dust surveys. The wavelength coverage of the spectro- 
scopic observations was tailored to include the [OI] 3 P 1 - 3 P 2 and 
3 P 0 - 3 P 1 lines at 63 and 145 pn t, [CII] 2 P 3 ^ 9 - 2 P 1 , 2 at 157 pm, 
several H 2 0 lines, particularly those at 78 and 180 ^m, along 
with adjacent transitions of CO and OH observable without in- 
curring a significant penalty on the total required time. In addi- 
tion, GASPS provides accurate far-infrared photometry at 70 
and 160 pm and, in most cases, 100 pm. M The project was 
awarded 400 hr of time to survey up to 250 young systems 
(in several cases, multiple systems were covered in the same 
observation), and observations were taken at various times be- 
tween 2010 December and 2012 July. The spectrometer was 
used with up to 4 settings per target, each of which covered 
a relatively small wavelength range (typically A A/ A ~ 5%) si- 
multaneously in two PACS grating orders. Most objects were 
observed in the setting covering the [OI] 63 pm line, with a sub- 
set of the brighter ones observed in the other settings, resulting 
in a two-phase survey strategy (see § 2.1.3). 

To help maintain the unbiased nature of the survey, targets 
were chosen with a wide range of spectral type, disk dust mass, 
age, and other stellar parameters (see § 3). They were located in 
seven well-studied young clusters and associations, with a dis- 
tance range of 40-200 pc (with the majority around 150 pc). 
Assuming typical disk sizes of 100-300 AU (Williams & Cieza 
2011) and with the angular resolution of Herschel/PACS of 5" at 
the shortest wavelength, line and continuum emission from the 
disk itself is unlikely to be spatially resolved. With the highest 
spectral resolution of PACS (88 km s _1 at 63 pm), disk emis- 
sion will also not be spectrally resolved (the Keplerian rotation 
velocity of most of the disk mass is ~10-50 km s _1 ). However, 
nondisk components such as outflow jets or ambient cloud 
emission may be resolved (see § 4). In most cases, all we have 
is a single measurement of the line flux on each target, yielding 
highly degenerate solutions to the underlying disk physics and 
chemistry. In the absence of resolution, a survey covering a wide 
range of target parameters is required, along with detailed 
modeling and data from other wavelengths. 

2.1. Observational Technique and Survey Strategy 

2.1.1. PACS Spectroscopy 

The PACS instrument (Poglitsch et al. 2010) offers resolu- 
tions of 1500-3400 (200-88 km s _1 ) and the ability to observe 
most of the 60-200 pm wavelength range. In spectroscopic 
mode, PACS provides an IFU with a 5 x 5 array of spaxels, 
and a pixel size of 9.2". By comparison, the instrumental point- 
spread function (PSF) ranges from 4.5" (FWHM) at 63 /im up to 
13" at 180 pm. For the GASPS project, spectral observations 


34 FIR fluxes for the brighter objects are available from IRAS, ISO or Spitzer, 
although in many cases the fluxes at A > 70 are unreliable because of the large 
beams and confusion levels involved — particularly in star forming regions. 


were taken using line-scan or range-scan modes, whereby the 
grating is scanned over a small wavelength range, taking data 
from all detector pixels simultaneously. Wavelengths around the 
central region of the spectra are observed by all of the 16 de- 
tector pixels to minimise flat-fielding problems due to inter- 
pixel variations. Line-scan observations have small wavelength 
coverage, and are designed to cover a single spectrally- 
unresolved line and immediately adjacent continuum with the 
full sensitivity. Range-scan observations have an arbitrary 
wavelength coverage and for GASPS were set up to include sev- 
eral close lines of interest by scanning up to 2 pm. Table 1 
shows the settings of the four wavelength scans A through 
D. Each has a primary line targetted in one of the grating orders; 
the secondary simultaneous grating order (given in brackets) 
was used to observe other useful lines (the full list of lines de- 
tected during the course of the survey is given in § 5.1). For a 
few individual targets, integration times longer that given in 
Table 1 were used for followup of marginal detections. Note 
that with the array spectral scanning technique, not all wave- 
lengths are being observed by a detector at all times, so the noise 
level increases towards the spectrum edges. However, the rms 
values in the scan centre in Table 1 are in good agreement with 
predictions. The observed sensitivity at 63 pm is equivalent to a 
3cr line luminosity sensitivity limit of 6 x 10 _6 (D/140) 2 L q , 
where D is the source distance in pc. 

Observations were performed in chop-nod mode with a small 
throw (1.5'), primarily to remove telescope and background var- 
iations. This chops out smooth background emission from 
around the source, but it may result in confusion from chopping 
onto extended emission. In some targets, this could be seen in 
the [CII] line (§ 4.5). 

2.1.2. PACS Photometry 

Photometric data were obtained using the fast scanning mode 
of the PACS imaging photometer, operating at central wave- 
lengths of 70 and 160 pm simultaneously (and repeated at 
100 and 160 pm in most objects). This technique scans the tele- 
scope over the source, using relatively short scan lengths of 3.5', 
and small (4") orthogonal steps between each scan. Two scans 
were performed, at 70° and 110° to the array, to improve the 
final image fidelity and avoid striping effects in the scan direc- 
tion. The photometer array field-of-view is 3.5' x 1.75' and, 
although the resultant image does not have a constant signal- 
to-noise ratio over the field, the noise level in the central 180" x 
80" region varied by less than 20%. The technique was found to 
be more sensitive than the chop-nod method, and in some cases 
several objects could be covered in the same field. It also en- 
abled searches for faint companions in the radius range 
~1500-7000 AU. The required sensitivity for the photometric 
observations was better than 5 mjy rms at 100 and 160 pm, and 
a factor of ~2 lower at 70 pm. Although the FIR continuum flux 
from disks is dependent on the stellar luminosity and mean disk 
temperature as well as the dust mass (and may be optically 
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TABLE 1 

PACS Wavelength Settings and Sensitivities in the Primary Grating Order 


Grating setting 

Primary line 

Primary A (/im) 

Grating orders 

Time 3 (s) 

Predicted rms b 

Observed rms b 

Notes 

A 

[OI]63.2 

63.08-63.29 

3 (1) 

1760 

4.6 

2.5-3 

LineScan 

B 

[CII] 157.7 

157.10-158.90 

1 (2) 

1500 

1.3 

0.8-1. 4 

RangeScan 

C 

H 2 0 180 

178.90-181.00 

1 (2) 

2000 

1.7 

2.0 

RangeScan 

D 

[OI] 145.5 

144.00-146.10 

1 (2) 

1630 

1.4 

1.6 

RangeScan 

Phot 

(Blue) 

60-85 

3 

180 

2.3 

2.6-3.0 


Phot 

(Green) 

85-130 

3 

180 

2.7 

2.6-3.0 


Phot 

(Red) 

130-210 

1 

360 

3.6 

4. 7-9.0 

(Background) 


a Approximate times for most targets in the main survey, not including overheads. 

b Note that noise rms levels are given in units of 10 -18 W m -2 for the spectroscopic observations, and in mJy for the continuum 
photometry. 


thick), disks of dust mass 10 -5 -10 -3 M G have IRAS 60 /im 
fluxes of typically 1 Jy at the fiducial distance of 150 pc. So 
the survey should detect dust in systems 1-2 orders of magni- 
tude fainter than this. The noise level of the observations was 
generally close to the original prediction (see Table 1), although 
in some cases it was limited by galactic background emission at 
160 /tm. Although the spectroscopic data could in principle be 
used to give narrowband continuum fluxes from the line-free 
parts of the spectra (albeit with a factor of ~30 less sensitivity 
than the full photometry), it was found that the photometric ac- 
curacy of these data was lower than the broad-band photometry, 
and generally they were not used for SED fitting. 

2.1.3. Phased Survey Strategy and Data Reduction 

The spectroscopic observations were carried out in one or 
more of the wavelength settings in Table 1. Phase I of the project 
consisted of [OI]63 /im observations of most targets using grat- 
ing setting A, concatenated observations of the brightest ~10% 
of targets in settings B-D (in order to reduce spacecraft slew 
overheads), plus photometry. 35 Note that not all objects from 
the initial survey list were observed in the lines; based on early 
survey results, a number of targets were dropped as they were 
deemed too faint in continuum to have likely emission in any 
line. In addition a few Taurus objects were dropped from both 
continuum and line observing based on updated re-classification 
as field stars (Luhman et al. 2009). Phase II of the project con- 
sisted of flexible followup of the brightest [OI] targets using 
grating settings B, C and/or D, as well as some deeper observa- 
tions of a few individual sources. 

During the course of the survey, GASPS photometric and 
spectroscopic data were reduced using prevailing versions of 
the standard Herschel data processing pipeline, HIPE (Ott 
2010). This provides calibrated FITS images and datacubes; 


35 Some bright targets were dropped from the photometric list, as suitable data 
was available from other Herschel surveys, e.g., some of the Taurus and Chall 
objects were covered by the Goult Belt Survey. 


further photometry and spectroscopic extraction were per- 
formed with packages such as STARLINK Gaia. However, the 
released version of HIPE evolved during the course of the mis- 
sion and different versions were used to reduce GASPS datasets 
in different publications, ranging from version 2.3 in early data 
(Meeus et al. 2010) to version 7.0 (Riviere-Marichalar et al. 
2012a) and version 9.0 (Howard et al. 2013, in preparation). 
Later HIPE releases generally have improved calibration as well 
as better flat-fielding, and the complete GASPS survey is to be 
re-reduced using a single mature version before being made 
publicly available as a systematically-calibrated dataset. The 
current work makes use of data extracted from the Herschel sci- 
ence archive during 2012, but the detection statistics presented 
here are unlikely to change significantly in the final data release. 

3. TARGET SELECTION 

GASPS targets were selected from the seven well-known 
nearby young star formation regions and associations listed 
in Table 2 and described in § 3.1. The complete target list, with 
system parameters from the literature, is given in Appendix A. 
For completeness, we list all the initial targets in this table, al- 
though some were not observed in spectral lines in the final sur- 
vey (see above). The criteria used to select the targets were: 

1. Age range 0.3-30 Myr. As discussed by several authors 
(e.g., Hartmann et al. 2001), stellar ages are uncertain — 
particularly for >10 Myr — and in these cases it may be better 
to take the ensemble age for a cluster rather than ages of indi- 
vidual stars. Systems of age <0.3 Myr were considered more 
likely to include nondisk emission components such as remnant 
ambient material, infalling envelopes, or energetic outflows (see 
§ 4). Those older than ~30 Myr were expected to have very 
little circumstellar gas. 

2. Optically visible stars. This means mostly SED Class II, 
III, Transition Objects, or debris disks. Targets have optical ex- 
tinctions less than ~3 m . We avoided embedded objects (i.e.. 
Class 0-1), because of potential confusion from extended sur- 
rounding gas. 
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TABLE 2 

Summary of Clusters and Associations in GASPS 


Group 

Distance (pc) 

Age (Myr) 

Disk fraction* (%) 

GASPS targets 

Notes/Main population 

Taurus 

140 

0.3^1 

90 

106 

Class I— III T Tauri stars 

Cha II 

178 

2-3 

75 

19 

Class II T Tauri stars 

77 Cha 

97 

5-9 

56 

17 

T Tauri and debris disks 

TW Hya 

30-70 

8-10 

>30 

13 

T Tauri and debris disks 

Upper Sco 

145 

5/11 

20 

44 

Class II— III T Tauri stars. 

t i Pic 

10-50 

10-20 

>37 

18 

Debris disks 

Tuc Hor 

20-60 

30 

>26 

16 

Debris disks 

HAeBe stars 

20-200 

-0.5-30 

100 

24 

Includes debris disks 


“Note: Disk fractions are based on published photometric excesses. They are the fraction of stars with any measured 
disk, so include mostly debris disks in the older associations and protoplanetary disks in the younger star forming 
regions. 


3. Disk dust masses 36 mostly in the range 10 -3 < M d < 
10“' M q . Also included were a number of coeval stars with 
Class III SEDs, or upper limits for M d of ~10 -6 M 0 , where 
the lack of continuum excesses in the IR or sub-mm suggested 
negligible warm or cool dust. Some of these still had gas accre- 
tion signatures and were included as they potentially could be 
associated with moderate masses of gas. In addition a number of 
debris disks in young associations were also observed, with M d 
as low as 10 -11 M 0 . 

4. Stellar spectral type AO through M 5. A similar range of 
stellar spectral types was chosen in each region where possible, 
although to increase the numbers of early spectral type objects, 
we also identified a sample of isolated well-studied Herbig 
AeBe stars with a similar age spread to that of the clusters. 
The resulting stellar mass range was ~0.2-3 M 0 , based on pub- 
lished HR diagrams. 

5. Nearby regions, with distances of <180 pc for low-mass 
and <200 pc for HAeBe stars. 

6. Low confusion level (from Herschel Confusion Noise 
Estimator). Confusion noise was <100 mly at 100 /tin. This 
meant that several dense star-forming regions such as p 
Ophiuchus were excluded from the survey. 

7. Extensive photometric and, in many cases, spectroscopic 
datasets available at other wavelengths. 

8. A range of accretion rates (based initially on Ha equiva- 
lent width, EW), X-ray luminosity, and binary separation. 

3.1. Individual Associations 

3.1.1. Taurus 

The Taurus star formation complex lies at 140 pc with a 
depth along the line of sight of ~20 pc (Torres et al. 2009). 
Taurus contains mostly low-mass stars with an age range of 
<0.1 up to ~10 Myr (Palla & Stahler 2002; Giidel et al. 2010). 


“Values of M d were mostly based on published nun wavelength continuum 
observations, with estimates based on shorter wavelength data in some cases. 


Stars in this region have been extensively studied at many 
wavelengths, and the census of Class 0-11 YSOs is essentially 
complete, with a significant fraction of the Class III YSOs also 
known (Rebull et al. 2010). Stellar parameters in Appendix A 
are from Kenyon & Hartmann (1995) and Kenyon et al. 
(1998), updated where appropriate by values in Furlan et al. 
(2006), Giidel et al. (2007) and Rebull et al. (2010), and with 
disk masses from Andrews & Williams (2005) and Currie & 
Sicilia- Aguilar (2011). 

3.1.2. Cha II 

Chamaeleon II is a nearby (178 pc; Whittet et al. 1997) star 
forming region included in the Spitzer Legacy cores-to-disks 
program (Evans et al. 2009) and the Gould Belt key program 
(Andre et al. 2010). It contains a lightly clustered distribution 
of low mass YSOs in a range of evolutionary states with spectral 
types K-M and ages estimated from protostellar evolutionary 
tracks of 4 ± 2 Myr (Spezzi et al. 2008). The GASPS subsam- 
ple consists of 19 targets from Cha II, generally with infrared 
colors of Class II objects and/or Ha equivalent widths indicative 
of gas accretion. Although there exists no deep sub-mm survey 
to give disk mass estimates (the survey of Young et al. [2005] 
only detected DK Cha and possibly IRAS 12500-765 8, with a 
dust mass limit on the other objects of 2 x 10 -4 M 0 ), Alcala 
et al. (2008) have estimated masses based on SED fits to far-IR 
data, and so we give these values (using their D01 models) in 
Appendix A. However, it should be noted that these entail higher 
uncertainties compared with estimates from the sub-mm. Spec- 
tral types are from Spezzi et al. (2008), binarity is discussed in 
Alcala et al. (2008), and L x comes from the ROSAT survey of 
Alcala et al. (2000). 

3.1.3. Upper Sco 

The mean cluster age and distance of Upper Sco was estimat- 
ed as 5 Myr and 145 pc by de Zeeuw et al. (1999), although 
more recent estimates suggest it may be as old as 11 Myr 
(Pecaut et al. 2012), which is consistent with the rather low disk 
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fraction observed. The 8-70 //m SED has been used to identify 
disks as Class II, Class III, or debris (Carpenter et al. 2009), and 
these classifications are given in Appendix A. The disk masses 
and system parameters are from the sub-mm observations of 
Mathews et al. (2012) and references therein. 

3.1.4. Tj Cha 

This is a compact grouping of ~19 stars, first identified as a 
young association through X-ray observations (Mamajek et al. 
1999). One of the reasons for interest in this cluster is its age, at 
an estimated 8 Myr, and relatively close distance (97 pc). The 
disk fraction, based on Spitzer observations at 24 and 70 /tm 
(Gautier et al. 2008), is 56%, which is relatively large for the 
age of the association (c.f., Haisch et al. 2001). A number of the 
stars are active accretors, and at least two are identified as Class 
II T Tauri stars (Sicilia-Aguilar et al. 2009). In the target list, 
spectral types, Ha EW, and the presence of an infrared excess 
are based on Sicilia-Aguilar et al. (2009), and information on 
binarity is from Bouwman et al. (2006). X-ray luminosities 
are taken from Mamajek et al. (1999) and Lopez-Santiago et al. 
(2010), and disk mass estimates are mostly based on FIR mea- 
surements (Currie & Sicilia-Aguilar 2011). 

3.1.5. TW Hya Association 

First recognised as a group of nearby young stars by Kastner 
et al. (1997), the number of members in the TW Hya association 
(TWA) is now at least 25 (Webb et al. 1999; Mamajek 2005). It 
is the closest association with accreting T Tauri stars, and in- 
cludes two classical T Tauri stars (TW Hya itself, and Hen 
3-600), and two bright debris disks (HD 98800, a hierarchical 
multiple system, and HR4796A). Low et al. (2005) used Spitzer 
to measure 24 and 70 /im excesses around TWA members and 
found these four systems have a 24 /im excess a factor of ~100 
larger than the other members. However, several of the other 
stars also have evidence of dust disks, from weak excesses at 
longer wavelengths. The age of this system is confirmed at 
~10 Myr (Barrado Y Navascues 2006). Parameters in Appen- 
dix A are taken from de la Reza & Pinzon (2004) and Scholz 
et al. (2007), with disk masses mostly from Matthews et al. 
(2007). Note that TWA member HR4796A (TWA11) is listed 
under the HAeBe stars as A-12. 

3.1.6. (3 Pic 

The moving group associated with /3 Pic was identified by 
Barrado y Navascues et al. (1999), and membership extended by 
Zuckerman et al. (2001) and others (see Torres et al. [2008] for 
a summary). With a derived mean age of 12 Myr, and range 
in distance of 10-50 pc, many of its members have been exten- 
sively studied over a wide range of wavelengths, including 24 
and 70 fim with Spitzer (Rebull et al. 2008), as well as the 
submm (Nilsson et al. 2009). The disk fraction is >37% (Rebull 
et al. 2008), and includes a number of debris disks in addition to 


/3 Pic itself. 37 Data in Appendix A are mostly from the above 
references. Detailed results from GASPS, including model fits 
to the photometry, have been presented for HD 181327, one 
of the brightest debris disks in this group, and HD 172555 
(Lebreton et al. 2012; Riviere-Marichalar et al. 2012b). 

3.1.7. TucHor 

This stellar association was first recognised by Zuckerman & 
Webb (2000) and Torres et al. (2000), who derived an age of 
20MB Myr and distance range of 20-60 pc. No N-band ex- 
cesses were seen around any stars in Tuc Hor (Mamajek et al. 
2004). However, a Spitzer study at 24 and 70 fim (Smith et al. 
2006) showed 5/21 stars with a measurable excess at 70 /im. 
Zuckerman et al. (2011) subsequently extended the search and 
found several more stars with IR excesses. All such systems in 
Tuc Hor are thought to be debris disks, and this is the oldest 
association in GASPS. The photometric data have been pre- 
sented in Donaldson et al. (2012). 

3.1.8. Herbig Ae Be stars 

The survey includes 25 IR-excess stars of spectral types late B 
to F, improving the statistics at the higher end of the stellar mass 
range (around 2-4 M 0 ). 38 This sample also includes some A-type 
stars with excesses where the classification is less clear: the pe- 
culiar Be star 51 Oph (Thi et al. 2005), and 5 systems which may 
be classified as debris or HAe (including 49 Cet, where the age 
was recently revised to 40 Myr [Zuckerman & Song 2012]). Like 
the lower-mass counterparts, the program HAeBe stars are biased 
toward isolated systems which have published IR excesses and 
ancilliary data (particularly UV spectra, resolved coronagraphic 
images and/or millimeter interferometry). HD 97048 — one of the 
brightest targets in our sample — had prior evidence from ISO of 
[OI] and possibly [CII] emission (Lorenzetti et al. 1 999). The HAe 
sample includes several disks with large gaps and/or cavities, as 
well as 2-3 systems with jets. Unlike the T Tauri stars, which are 
represented in sufficient number to permit statistical evaluation of 
association ages, the HAe stars represent extremes in stellar and 
disk properties, and have more uncertain ages except where there 
are common proper motion late-type companions. Stellar para- 
meters in Appendix A are mostly taken from Montesinos et al. 
(2009) and Meeus et al. (2012), with disk masses from Acke et al. 
(2004) and Sandell et al. (2011). Note that the SED classifications 
in the appendix are different from the T Tauri class, and are based 
on the mid-IR slope as suggested by Meeus et al. (2001): Group I 
has an SED rising to longer wavelengths in the mid-IR, and group 
II has a falling SED. Results from the GASPS HAeBe subsample 
have been presented in Meeus et al. (2012). 


37 /3 Pic was observed as part of the Herschel GT Program ‘Stellar Disk Evo- 
lution' (P.I. G. Olofsson). 

38 Note that the HAe star AB Aur is listed under the Taurus subsample as T-101 
in Appendix A. 
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3.2. Ancilliary Data 

Many of the GASPS targets are well-known systems, with 
photometry in optical through to mid-IR (including Spitzer 
fluxes at wavelengths as long as 70 fi m), and sub-mm (mostly 
850 / /m). In addition, Ha or Br 7 line strengths are published for 
many targets. Derived parameters such as stellar spectral type, 
T eff , disk dust mass, and SED Class are also mostly available, 
and the most recent published estimates are given in Appen- 
dix A. As part of the GASPS project, we have endeavoured to 
obtain such data in cases where it is missing, and to make the 
target sample uniform both in data and in derived parameters. 
One additional issue is that much of the published photometry is 
not contemporaneous; in some cases, photometric points in the 
optical and NIR have been taken 20 years apart. For time- 
variable objects, SED fitting under these circumstances may 
be significantly affected, and more recent optical and near-IR 
photometry is being obtained for a number of the targets in order 
to improve the reliability of SED fits. 

4. ORIGINS OF FAR INFRARED LINES FROM 
YOUNG STARS 

At far-infrared wavelengths, common species such as C, O 
and N have several prominent fine-structure transitions. [01] 
lines at 63 and 145 /im and the [CII] line at 157 /im are impor- 
tant cloud coolants, on a galactic scale (Stacey et al. 1991), in 
photodissociation regions (Hollenbach et al. 1991), and in cir- 
cumstellar disks (Kamp et al. 2003; Gorti & Hollenbach 2004). 
In star formation regions, both models and observations indicate 
that [OI]63 /rm is the single brightest emission line in the FIR/ 
sub-mm. Abundant molecules such as CO and H 2 0 also have 
numerous rotational lines throughout the FIR with energy levels 
of a few hundred K, and can trace the ‘warm’ gas components. 
Other FIR-emitting species such as OH are photodissociation 
products of H 2 0, and are therefore predicted to be abundant 
(Najita et al. 2010). Around individual young stars, FIR lines 
can arise from several different regions. For example, CO and 
OH emission from young highly luminous HAeBe stars was 
thought to arise from dense regions of size ~200 AU (Giannini 
et al. 1999), but it was unclear whether these were disks or rem- 
nant envelopes. Lines are also seen from high-velocity jets and 
low-velocity photoevaporating disk winds, and the relative con- 
tributions of disk, outflow, disk wind, and envelope will depend 
on the SED class, stellar radiation field, disk structure, mass loss 
rate and the environment. 

The GASPS project involves both in-depth studies of indi- 
vidual targets using multi-wavelength data, as well as a statisti- 
cal analysis of the full FIR sample. Interpretation of the results 
generally requires detailed comparison with models, and in the 
following sections we outline methods of estimating the contri- 
butions to the FIR line emission, focussing on Class II— III 
YSOs, which form the bulk of our targets. 


4.1. Disks 

Fine structure atomic line emission arises from the surface of 
disks at Ay ~ 1 over a wide range of radii, where the stellar UV 
or X-ray photons ionise the exposed gas to produce a mainly 
atomic extended disk atmosphere (Meijerink et al. 2008; Gorti 
& Hollenbach 2008; Woitke et al. 2009a). [OI]63 //in is pre- 
dicted to be the brightest line from disks at any wavelength, with 
line luminosities as high as ~10 -4 L 0 from T Tauri systems 
(Gorti & Hollenbach. 2008). It becomes optically thick relative- 
ly easily, and traces the mean gas temperature on the disk sur- 
face rather than the mass directly. FIR molecular lines such as 
CO and H 2 0 also arise from the warm heated surface of dense 
disks (Woitke et al. 2009b). However, molecular photodissoci- 
ation in more tenuous debris disk systems may mean the atomic 
fine structure lines will dominate the FIR (Kamp et al. 2003; 
Zagorovsky et al. 2010). Line fluxes depend strongly on the disk 
structure (for example, a flared disk has a larger exposed surface 
area, resulting in brighter lines [Jonkheid et al. 2004]), the ra- 
diation field from the central star, as well as the details of chem- 
istry and gas/dust ratio in the disk atmosphere. In the following 
we summarise the methods used for modeling the data. 

4.1.1. Disk Modeling 

MCFOST and ProDiMo are the two main codes used in 
GASPS to model the protoplanetary disk structure and appear- 
ance. MCFOST is a three-dimensional Monte Carlo continuum 
and line radiative transfer code (Pinte et al. 2006, 2009). The 
parametrized input disk density distributions can accommodate 
structures such as holes, gaps and dust settling. The calculation 
of the dust temperature and radiation field takes into account 
non-isotropic scattering, absorption and re-emission based on 
the local dust properties. The code uses a large variety of grain 
size distributions and compositions, e.g., porous grains and icy 
grains. SEDs, thermal and scattered light images, visibilities, as 
well as line emission are derived by a Monte-Carlo method and 
ray-tracing of the final physical disk structure. ProDiMo is a 
two-dimensional thermo-chemical disk code that calculates the 
vertical hydrostatic equilibrium, gas phase (e.g., neutral-neutral, 
ion-molecule, photochemistry. X-ray chemistry), and gas-grain 
chemistry (ad- and desorption processes), two-dimensional con- 
tinuum radiative transfer (with isotropic scattering), detailed gas 
heating/cooling processes (including X-rays) using 2-directional 
escape probability, and spatial decoupling of gas and dust (e.g., 
settling) (Woitke et al. 2009a; Kamp et al. 2010; Aresu et al. 
2011; Woitke et al. 2011). The observables derived from the 
resulting chemo-physical disk structure include SEDs, line 
fluxes, profiles, and images. For optically thin cases such as de- 
bris disks, we also use GRaTer, a ray-tracing code incorporat- 
ing a large variety of grain compositions which fits SEDs, 
images, and interferometric visibilities using parametrized op- 
tically thin disk models (Augereau et al. 1999). 
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r [AU] r [AU] 


Fig. 1. — Cross-sections through ProDiMo models of a T Tauri disk (left) and a HAeBe disk (right), with the density structure as a greyscale, and the Ay = 10 m 
surface shown as a dashed red line. The primary emitting regions for GASPS lines are given by the coloured boxes, and indicate where 80% of the emission arises. The 
lines are [OI]63 /nn (green), CO J = 18 — 17 (orange), [CI1J157 //rn (black), OH79 //m (blue) and the 63.3 /mi H 2 0 line (white). Also shown is the CO J = 3 — 2 
emitting region (red). Model parameters are given in Table 3. Note the larger outer radius displayed in the HAeBe model. 


For modeling individual sources, grids of MCFOST or 
GRaTer models were run over a broad parameter space to iden- 
tify the best fitting dust model, based on SEDs, images, and 
interferometric data when available. MCFOST results were 
passed to ProDiMo for detailed gas modeling. Examples of this 
approach are Meeus et al. (2010) and Thi et al. (2010) and for a 
debris disk, Lebreton et al. (2012). Another approach employs a 
genetic algorithm minimisation strategy with the ProDiMo 
models to find local minima in the parameter space constrained 
by observations. Examples are found in Tilling et al. (2012) and 
Woitke et al. (2011). 

In Figure 1 we use ProDiMo to illustrate the regions where 
most FIR line emission is expected to arise in T Tauri and 
HAeBe disks, using the model parameters given in Table 3. 
Note that these are relatively massive disks — towards the high 
end of the range of the GASPS sample. For these models we 
have proscribed the disk vertical struture by the scale height 
and flaring index. The results indicate that the [OI] lines are seen 
mostly from the disk surface at 20-200 AU radius around a 
T Tauri star, and a factor of 1 .5 further out in the more luminous 
HAeBe star. The [CII] line comes from the tenuous outer atmo- 
sphere at radii >100 AU, whereas high-J CO emission (e.g., the 
J = 18 — 17 line at 144.78 /im) is predicted to arise only from 
within a few tens of AU for T Tauri disks. The model predicts 
emission from the 63.3 urn line of H ,0 mostly from within a 
few AU. 

4.1.2. Model Grids 

To support a broader statistical analysis of the GASPS data, 
we have produced a grid of models covering a wide parameter 


space (Woitke et al. 2010). 39 Stellar masses between 0.5 and 
2.5 M e and pre-main sequence evolutionary tracks at 1, 3, 
10, and 100 Myr are used to define T eff , H t and hence L t . 
The UV excess, / uv , is treated as a power-law that is added 
on top of the photospheric spectrum in the wavelength range 
912-2500 A. Dust masses range from 10 -7 to 10 -3 M 0 , and 
the gas/dust mass ratio runs fromlO 3 (10 x ISM) to 0.1 (0.001 x 
ISM). Geometries include young flaring disks, flat evolved sys- 
tems, as well as ‘transition’ disks with inner holes up to 100 
times the sublimation radius. The grid also contains models with 
a settled dust distribution, where larger grains have a smaller 
vertical scale height than smaller ones. The observables calcu- 
lated from these models include SEDs and integrated line 
fluxes. 

Woitke et al. (2010) show that the fine structure line fluxes of 
[OI] and [CII] depend strongly on the stellar UV excess and disk 
flaring. Using the entire grid of parameter space (not folding in 
the likelihood of these disks occurring in nature), about 70% 
of the high-mass models (dust mass, M d > 10 -5 M 0 ) with a 
strong UV excess were predicted to be detected in [OI]63 //in 
line by GASPS, and 55% detected in [CII] 157 /im. Without a 
UV excess, the percentage drops to ~30% (14% for [CII]). An 
initial statistical comparison between the early GASPS line 
fluxes and the grid (Pinte et al. 2010) shows that some of the 
disks around low-mass stars (<2 L 0 ) do require additional 
heating from a moderate UV excess (with / uv = 0.1) or X-rays 
(which were not included in this first grid). However, results 


39 The grid was calculated on the FOSTINO computer cluster financed by 
ANR and operated by SCCI at OSUG. 
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TABLE 3 


Parameters Used by ProDiMo to Illustrate the Region of Line 
Emission from Disks Around T Tauri and HAeBe Stars in Fig. 1 


Parameter 

T Tauri 

HAeBe 

Notes 

SpT 

K4 

A3 


M, (M 0 ) 

1.4 

2.2 


/ 'K:- 

4500 

8600 


L* (L Q ) 

2.0 

32.0 


T (Myr) 

2.3 

4.6 

Age 

/uv 

0.01 

0 

Additional UV fraction, 
L\j\/ T* 

Lx (erg/s) 

10 3n 

0 

Additional X-ray 
luminosity 

Rin (AU) 

0.1 

0.5 

Set by the grain 
sublimation radius 

Rout (AU) 

300 

500 

Outer disk radius 

M d (M 0 ) 

10~ 4 

io- 4 

Disk dust mass 

g/d 

100 

100 

Gas/dust mass ratio 

e d 

1.0 

1.0 

Surface density power 
law exponent 

^min/^max (/^^n) ■ ■ • • 

. 0.05/1000 

0.05/1000 

Min/max grain size 

p 

3.5 

3.5 

Grain size power law 

/ PAH 

0.01 

0.01 

PAH mass fraction 

H„ (AU) 

10 

10 

Scale height at 100 AU 
radius 

0 

1.1 

1.1 

Disk flaring index 


from T Tauri disk modeling with stellar X-rays indicate that the 
[OI]63 /im line flux is only affected by X-ray heating for L x > 
10 30 erg/s (Aresu et al. 2011). 

Although most gas emission lines sample a thin, warm sur- 
face layer (see Fig. 1), combining the FIR data with results from 
other wavelengths (e.g., the [OI]63 /im/CO(2-l) ratio) and with 
physically-plausible models does allow us to break model de- 
generacies, giving approximate estimates for gas parameters in- 
dependent of the dust (Kamp et al. 201 1). But it is clear that the 
reliability of derived values such as the gas mass relies on the 
accuracy of the models. 

4.2. Outflow Jets 

Highly embedded Class 0-1 YSOs are known to have prom- 
inent outflows, and early observations using the KAO as well as 
more recent observations with Herschel/PACS show bright FIR 
lines around Herbig-Haro objects and high-velocity CO outflow 
lobes (Cohen et al. 1988; van Kempen et al. 2010). As well as 
fine structure lines such as [01] 63 /mi, many CO and H 2 0 tran- 
sitions are readily detectable from Class 0-1 objects (Lorenzetti 
et al. 2000; Molinari et al. 2000; van Dishoeck 2004; 
Goicoechea et al. 2012). The stars are young (<0.1 Myr), op- 
tically obscured, embedded in an envelope, and located near 
dense cloud cores. Their dense environments and high outflow 
luminosities suggest that the FIR line emission is dominated by 
outflow shocks (Molinari et al. 2000; Nisini et al. 2000, 2002; 
Franklin et al. 2008). These shocks also affect molecular 


abundances; for example, releasing H 2 0 from grains and in- 
creasing its gas-phase abundance to as much as 10 -4 — 
comparable with that of CO (Benedettini et al. 2002). In most 
cases the FIR lines dominate the shock cooling, and line fluxes 
may be used to estimate the outflow luminosity (Hollenbach 
et al. 1985; Nisini et al. 2002; Podio et al. 2012). 

Evolved, isolated objects such as optically-visible Class II- 
III T Tauri stars have mass accretion rates at least 1-2 orders of 
magnitude lower than Class I objects (Hartmann et al. 1998; 
Arce & Sargent 2006), which are themselves an order of mag- 
nitude lower than the Class 0 objects (Bontemps et al. 1996; 
Podio et al. 2012). Class II objects are pure disk systems, and 
are generally isolated with no dense ambient gas. Consequently 
we assume the fraction of outflow luminosity deposited in 
shocks near the star, r/ s , is given by the geometric fraction of 
the initially broad wind intercepted by the disk. So r/ s ~ H 0 / 
100, where H 0 is the scale height (in AU) at 100 AU radius. 
In the same way as embedded objects, the jet mass loss 
rate M out can be estimated from the [OI]63 /im luminosity, 
L 0 1 by: 

Mout \Vw'f 0\'Vs\^a/ ^w)\ i 

where f ol is the fraction of FIR line luminosity in the 
[OI]63 /tm line, v w the outflow jet velocity, and v a the ambient 
shock velocity (Nisini et al. 2002). For embedded objects, the 
dominant emission is from the integrated CO and H 2 0 lines, 
thought to be from C shocks, and Goicoechea et al. (2012) 
find /oi ~ 0.12 in the Class 0 YSO Serpens SMM1. In fast dis- 
sociative J shocks, [OI]63 /im emission may dominate the lu- 
minosity, and /oi is found to be 0.5 or greater (Podio et al. 
2012). Assuming a jet velocity of 100 kms -1 (Podio et al. 
2012), v a /v w = 0.2 (Nisini et al. 2002), with H 0 = 10 AU 
(Table 3), then the GASPS sensitivity (§ 2.1.1) may allow the 
detection of outflow mass loss rates of ~3 x 10 -9 M Q /yr for 
stars at a distance of 140 pc. However, if C shocks dominate, 
the contribution to the [OI] line from the jet may be lower. In 
Appendix A we have indicated the stars with published mass loss 
rates greater than this value (Hartigan et al. 1995). 

Although most are isolated disks, a few of the GASPS targets 
are somewhat embedded Class II objects and have extended op- 
tical jets; for these we may expect some outflow contribution to 
the FIR lines. The spatial resolution of PACS is modest, but can 
help investigate this contribution, resolving jets on scales of 
>1000 AU. For spatially-unresolved outflows the situation is 
less clear. However, shocked outflow emission may be broad- 
ened to ~200 km s _1 or velocity-shifted by more than a few 
tens of km s _1 , similar to the high-velocity component seen 
in optical lines in a few high-accretion objects (Hartigan et al. 
1995; Acke et al. 2005). In these cases, the PACS spectral reso- 
lution of 88 km s _1 at 63 /tm may also be used to help discrimi- 
nate between outflow and disk. 
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4.3. Remnant Envelope Gas 

Low-density PDRs in the remnant envelope gas centred on 
the stars may contribute to the [CII] flux from some objects, as 
the [CII] 157 pm critical density is only ~3 x 10 3 cm l To mit- 
igate this, targets were selected to be SED Class II — III with low 
A v , i.e., optically-visible stars where the envelope mass is at 
least 1 or 2 orders of magnitude less than the disk mass (Fuente 
et al. 2002; Arce & Sargent 2006). The [OI]63 pm critical den- 
sity is ~ 100 times higher than [CII], and the mean envelope den- 
sity is small compared to the disk, so the envelope contribution 
to the total [OI] flux should be small. 

4.4. Disk Winds 

A photoevaporative UV-driven wind (Pascucci & Sterzik 
2009) will serve to extend the scale height of a disk atmosphere, 
and may enhance FIR emission lines. The effect of this on the 
[OI]63 pm flux is under investigation. 

4.5. Extended Ambient Gas 

Observations of star-forming clouds in FIR fine structure 
lines using the LWS on ISO showed bright emission in regions 
containing luminous Herbig AeBe and FU Ori stars (Lorenzetti 
et al. 1999, 2000; Creech-Eakman et al. 2002). The [CII] emis- 
sion in many of these objects was spatially extended, and com- 
monly the line fluxes at reference positions many arcminutes 
from the stars were as bright as the on-source position. In these 
cases, the dominant [CII] source was thought to be extended 
low-density PDRs, with optical depths of l-2 m (Hollenbach 
& Tielens 1997). The GASPS targets were chosen to avoid 
the densest clouds, and we used the Herschel Confusion Noise 
Estimator (HCNE) to estimate the 100 pm continuum confusion 
noise (F c ) for all targets. From this we estimate a [CII] confu- 
sion noise level by adopting a ratio of T[cn|/-U c = 1.2x 
10 -19 Wm _2 /mJy, found for large-scale Galactic clouds by 
Shibai et al. (1991). During the course of the survey, 10 objects 
having a relatively high continuum confusion level ( F c > 
30 mjy) were observed in [CII] (8 of which were in the Taurus 
cloud). Based on the HNCE, the predicted [CII] confusion level 
for these was >3.6 x 10 _ls W m~ 2 . An examination of the ini- 
tial data shows no extended [CII] over the PACS footprint in 9 of 
these objects, with an rms level of ~2 x 10 -18 Wm -2 . Either 
the confusion level is lower than predictions from the HCNE, 
or the [CII] emission is smooth over the PACS IFU field 
(~arcmin) and emission is being chopped out. One high back- 
ground confusion source (HD 163296) had evidence of extend- 
ed [CII] at a level of ~10 -17 W m -2 in the PACS field of view 
and in the chopped reference beam. From the HCNE, this object 
has the highest value of F c in the GASPS sample (85 mJy), 
which would predict, based on the above ratio, a [CII] confusion 
noise of J[ CII ] = 10“ 1 ' Wm' 2 , consistent with the observations. 

The lOOx larger critical density of [OI]63 pm compared with 
[CII] implies that extended [OI] emission from diffuse ambient 


gas is expected to be negligible (Liseau et al. 2006). ISO found 
that the 63 pm line flux is mostly higher towards highly lumi- 
nous YSOs than off-source. 

4.5.1. Line-of-Sight Absorption 

As well as emission, dense clouds may have significant op- 
tical depth and be self-absorbed in the [OI]63 line. However, 
estimates suggest the line optical depth may not become signif- 
icant until Av > 10 (Liseau et al. 2006; Abel et al. 2007). 
Moreover, the linewidths of the cool line-of-sight clouds are 
<1 km s _1 , small by comparison with the 5-20 km s _1 widths 
predicted for Keplerian disk emission. Combined with the ex- 
tinction limit of Av < 3 in the GASPS survey means this effect 
should be small in most cases. 


5. FIRST RESULTS 

Results from some subsets of the GASPS study have been 
presented in previous papers. A summary of the ‘science dem- 
onstration' observations of a small number of targets was given 
in Mathews et al. (2010), and a comparison of these data with a 
broad grid of disk models was shown by Pinte et al. (2010). 
More detailed comparisons of the line and continuum data with 
individual tailored models were carried out based on the detec- 
tions of [OI]63 pm in the T Tauri star TW Hya (Thi et al. 2010) 
and the HAeBe stars HD 169142 (Meeus et al. 2010) and 
HD 163296 (Tilling et al. 2012) . The T Tauri star ET Cha 
was detected in both [OI]63 pm and FIR continuum, and model- 
ing indicates the disk is unusually compact (Woitke et al. 2011). 
CH + was detected in one of the brightest targets — the HAeBe 
system HD 100546 (Thi et al. 2011). An emission line of H 2 0 at 
63.3 pm found in a number of the T Tauri stars indicates warm 
(~500 K) H 2 0, possibly from the inner few AU of the disks 
(Riviere-Marichalar et al. 2012a)). In most of the older gas-poor 
systems the lines were not detected; however, the far-IR pho- 
tometry has been used to improve the SEDs and dust modeling 
(Donaldson et al. 2012; Lebreton et al. 2012). 

In the following sections we summarise some of the overall 
results from GASPS, including identification of the lines found 
in the survey, and an initial comparison of the spectra of different 
types of objects (§ 5.1). In Appendix A (column 12) we indicate 
which of the four primary lines ([OI]63 pm, [CII] 157 pm, 
CO J = 18 — 17, and H 2 0) were observed and detected 
in the targets. For these purposes, a detection is regarded as 
>3-<r above the noise. In § 5.2 we give the overall line detection 
statistics from the survey, and discuss the [OI]63 pm and [CII] 
emission characteristics in § 5.3 and 5.4. Finally, in § 5.5, we 
show the effects of other system parameters on the line detection 
statistics. It should be noted that these data are mostly based on 
results from early versions of the reduction pipeline HIPE. Con- 
sequently, the flux calibration and flat-fielding is not finalised 
and some detections are subject to re-analysis. Final values of 


2013 PASP, 125:477-505 


This content downloaded from 128.154.208.61 on Tue, 9 Jul 2013 16:21:27 PM 
All use subject to JSTOR Tenns and Conditions 


488 DENT ET AL. 


the fully-calibrated fluxes and detailed flux correlations will be 
given in subsequent papers. 

5.1. Summary of Lines Detected 

To illustrate and compare the lines detected in the richer 
GASPS targets, spectra from the central spaxel in three objects 
from the survey are shown in Figures 2-5 (note that the spectra 
are scaled to enable comparison in these figures). T Tau (Podio 
et al. 2012) ( shown in red) is a KOV star with a massive disk, 
compact outflow, some surrounding reflection nebula, and pos- 
sibly a PDR. FIR lines may arise from a mixture of these com- 
ponents. although the molecular transitions seen in the ISO LWS 
spectra were attributed mainly to shock emission (Spinoglio 
et al. 2000). HD 100546 (Meeus et al. 2012) ( shown in blue) 
is a young Herbig AeBe star with a bright disk but without a 
prominent outflow, but which also has a rich FIR spectrum. 
AA Tau (in green) is perhaps a more typical isolated T Tauri 
star, with a luminosity of ~1 L 0 , weak outflow and a relatively 
massive disk. In Table 4 we identify all the lines observed in 
these three objects. 

Both T Tau and HD 100546 have similar strengths in the 
fine-structure atomic lines. AA Tau is ~200 times weaker, but 
is detected in [OI] with a comparable line/continuum ratio to the 
others. However, it shows no evidence of [CII] . In HD 100546, 




microns 

Fig. 2. — Spectra of HD 100546 (blue), T Tau (red) and AATau (green), taken 
from the central PACS spaxel in the two shortest wavelength observations. 
Fluxes of T Tau and AA Tau fluxes are scaled by 2 and 150 respectively to 
facilitate comparison of the spectra. Lines found in any of the datasets are iden- 
tified (although not all the lines are seen in all objects) — see Table 4 for full 
details of the transitions. Note that small wavelength errors are sometimes ap- 
parent in these early reductions of the AA Tau spectra. 




Fig. 3. — Same as Fig. 2, for regions around 79 and 90 /in l. The spectra of 
T Tau and AA Tau have been multiplied by 2 and 150. The emission close to 
90 fim is a blend of H 2 0 and CH + , and in HD 100546 is thought to be mostly 
fromCH+. 


molecular transitions have a line/continuum ratio which is con- 
siderably lower than both T Tau and A A Tau. The PACS data 
cover four transitions of CO: J = 18 — 17, 29-28, 33-32 and 




Fig. 4. — Same as Fig. 2 for regions around 145 and 158 fim. The spectra of 
T Tau and AA Tau have been multiplied by 2 and 150. 


2013 PASP, 125:477-505 


This content downloaded from 128.154.208.61 on Tue, 9 Jul 2013 16:21:27 PM 
All use subject to JSTOR Terms and Conditions 


HERSCHEL GASPS— SUMMARY AND STATISTICS 489 



Fig. 5. — Same as Fig. 2 for region around 180 //in. The spectra of T Tau and 
AA Tau have been multiplied by 2 and 150. 


36-35. T Tau shows emission in all four CO lines, and compar- 
ison with the CO rotational diagram of Spinoglio et al. (2000) 
shows that the three highest transitions are new detections, re- 
quiring an additional hot gas component (>1000 K) to account 
for the emission. AATau is detected only in the two lower-level 
CO lines, most likely because of sensitivity limits. 

The OH doublet around 79 /an is detected in all three sources 
(Fig. 3). Several H 2 0 lines with upper energy levels from 1 15 to 
1300 K are seen towards both T Tau and AA Tau, and in AA 
Tau, H 2 0 is the only line detected, other than [OI]63 /tm, CO 
and OH. By contrast, HD100546 has no clear evidence of H 2 0 
emission, although other lines (atomic species, OH and C H ) 
are relatively bright. The highest energy level H z O transition 
covered by GASPS is the ortho 8 18 -7 0 7 line at 63.3 /an (Fig. 2); 
this was detected in T Tau, AA Tau and several other T Tauri 


stars (Riviere-Marichalar et al. 2012a)). Finally, both HD 
100546 and T Tauri show clear CH + emission at 72.14 /im, with 
possible blends of CH + and H 2 0 around 90.0 and 179.5 //m; 
this species was also identified at several other wavelengths in 
HD 100546 (Thi et al. 2011). 

5.2. Primary Line Detection Statistics 

Column 12 of Appendix A shows whether each of the four 
primary species ([OI]63 /an, [CIIJ157 /an, CO J = 18 — 17 
and H 2 0 63.3 //in) were detected in the GASPS targets. Based 
on this, the overall detection statistics are given in Table 5. As 
noted above, observations of the range-scan observations were 
normally only performed if a target was already found or ex- 
pected to be detected in [OI] 63 /m i, so the detection rates of 
[OI] 145 /xm, [CD], and CO in this table are biased towards those 
with known [OI]63 /tm emission. Of targets observed in multi- 
ple lines, only one remained undetected in [OI]63 /an yet shows 
emission in one of the other lines. Based on this result and our 
modeling, it is thought unlikely that a significant number of the 
[OI] -unobserved objects would show emission in these other 
lines. The H 2 0 rates are the fraction of targets seen at 63.3 //in, 
which was observed as part of the [OI]63 /xm line-scan obser- 
vations. The CO rates are the fraction of targets detected in the 
brightest line covered by GASPS (CO J = 18 — 17). 

The main similarities and differences between line emission 
from the two types of objects are: 

1. Of the sample of 164 objects observed in spectroscopy at 
63 /xm, approximately 49% were detected in [OI]. 


TABLE 4 


Lines Identified in HD 100546, T Tau and/or AA Tau 


Wavelength (/rm) 

Line ID 

Transition 

^Tipper (K) 

HD 100546 

T Tau 

AA Tau 

63.18 

[OI] 

3P1-3P2 

228 

X 

X 

X 

63.33 

o-H 2 0 

$18 — 7()7 

1293 


X 

X 

71.94 

o-H 2 0 

^07 6 16 

685 


X 

X 

72.14 

CH+ 

7 = 5-4 

600 

X 



72.84 

CO 

J = 36 - 35 

3700 


X 


78.74 

o-H 2 0 

423 — 3i2 

432 


X 

X 

78.92 

p-h 2 0 

@15 — ^24 

396 


X 


79.11 

OH 

1/2 - 3/2 hfs 

182 

X 

X 

X 

79.36 

CO 

7 = 33 - 32 

3092 

X 

X 


89.99 

p-h 2 0 

322 — 2n 

297 

(Blend with CH + ) 

X 


90.02 

CH+ 

7 = 4-3 

297 

X 



90.16 

CO 

7 = 29 - 28 

2400 

X 

X 

X 

144.52 

p-h 2 o 

4l3 — 322 

396 


X 


144.78 

CO 

7 = 18-17 

945 

X 

X 

X 

145.52 

[OI] 

3PO-3P1 

326 

X 

X 


157.74 

[CII] 

2P3/2 - 2P1/2 

91 

X 

X 


158.31 

p-h 2 0 

331 — 4o4 

410 




179.53 

o-H 2 0 

2l2 ~ lot 

115 

(Blend with CH + ) 

X 


179.6 

CH+ 

7 = 2-1 

114 

X 



180.49 

o-H 2 0 

221 — 2i2 

194 


X 



An A' indicates a detection. 
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TABLE 5 

Detection Statistics of Primary Atomic and Molecular Species 
[OI]63 [OI]145 [CII]157 H 2 0 63 CO 18-17 


Total 80/164 24/61 19/72 12/164 24/58 

HAeBe stars* 20/25 5/23 6/25 2/25 10/24 

T Tauri stars' 1 60/139 19/38 13/47 10/139 14/34 


Each entry gives the number of targets detected and number observed. For 
[OI] 145, [CII], and CO, observations were mostly carried out only if the lines 
were detected (or likely to be detected) in [OI]63 /im. 

“ Includes five young A stars classed as debris disks. 

b This includes all stars observed which were not part of the HAeBe group. 


2. A biased subset of the brighter objects from (1) were ob- 
served in [OI] 145 /im, [CII] 157 /im and CO J = 18 — 17, and 
the detection rates in this subset were 25^10% in each of these 
lines. Assuming that [OI]63 /im is always the easiest to detect 
(see above), then an unbiased sample of all 164 targets from (1) 
would have had a detection rate of ~14% in these other lines. 

3. All HAeBe stars were detected in [OI]63 pm — a signifi- 
cantly higher detection rate than T Tauri systems. (Note that the 
statistics of HAeBe stars in Table 5 include five known A-star 
debris disks.) 

4. The [OI]145 detection rate is a factor of ~2 higher in the 
T Tauri stars observed compared with HAeBe systems. This 
may reflect a higher [OI]63/145 pm line ratio in HAeBe disks. 

5. The [CII] detection rate is similar (26%) in both T Tauri 
and HAeBe stars. If this is envelope material (see § 5.4), it in- 
dicates that compact envelopes of atomic gas can be maintained 
around both high and low-luminosity stars. Note, however, 
that in some cases the [CII] emission may be confused by am- 
bient gas. 

6. One (possibly two) HAeBe stars were detected in H 2 0. 
Although in the small number regime, the H 2 0 detection rate 
is formally similar to that of T Tauri systems. However, consid- 
ering the HAeBes are relatively bright in continuum compared 
with the T Tauri sample, this suggests that, on average, HAeBe 
systems are weaker in H 2 0 compared with T Tauri systems. 

7. The fraction of objects with detectable warm CO (based on 
the J = 18 — 17 transition) is similar (40%) in disks around 
both types of stars. 

5.3. [OI] Line Emission 

As is clear from the example spectra, [OI]63 pm is normally 
several times brighter than any of the other FIR lines observed 
by GASPS, with an overall detection rate in the survey of ~49%. 
In most cases, it is the best tracer (in the far-IR) of whether gas is 
present. This is true for almost all GASPS sources. To help un- 
derstand the origin of the emission we can look at the data in 
more detail. 

Most objects were unresolved in both line and continuum 
emission. An example is AA Tau (Fig. 6), where the ratio of 
flux in the centre to average of adjacent spaxels is ~20. This 



Fig. 6. — Spectra covering the [OI]63 /tm line from the compact source AA 
Tau, in the central spaxel ( red histogram ), and an average of the eight adjacent 
spaxels (in blue). The adjacent pixel spectrum has been scaled up by a factor of 5 
for clarity. Both the [OI]63 //in and nearby I I 2 0 line are detected only in the 
central spaxel. 


is consistent with an unresolved source, where we would expect 
the adjacent pixel average to be a few % of the centre, given an 
inter-spaxel spacing of 9.4", a PSF Gaussian equivalent width of 
~5.4" at 63 /tm, the asymmetric sidelobes from PACS of a few 
percent, and taking into account possible pointing uncertainties 
of a few arcsec in some datasets (PACS User Manual, 2011). 
This lack of extended emission indicates a line emitting region 
of radius <500 AU. 

For a number of individual unresolved objects with low ac- 
cretion rates and no evidence of outflow we have assumed a disk 
origin and combined the [OI] fluxes with data at other wave- 
lengths to estimate disk properties. Initial ProDiMo models 
of the relatively large disk in TW Hya (several 100 AU radius) 
indicate a gas mass of a few 10 -3 M 0 with gas-to-dust ratio a 
factor of ~10 lower than the ISM value (Thi et al. 2010), al- 
though some models suggest the gas mass an order of magni- 
tude larger, with a more ISM-like gas-to-dust ratio (Gorti et al. 
2011). ET Cha, by contrast, has a compact disk of modeled ra- 
dius of only 10 AU, a low dust mass of a few 10 -8 M 0 and gas 
mass of a few 10~ 4 M 0 (Woitke et al. 2011), suggesting either 
the gas-to-dust ratio is enhanced or there may be another con- 
tribution to the line flux. The HAeBe stars HD 169142 and 
HD 163296 both show emission consistent with disks and 
ISM-like values of the gas/dust ratio (Meeus et al. 2010; Tilling 
et al. 2012). 

5.3.1. Spatially and Spectrally Resolved [OIJ63 pm 
Emission: Outflow Jets 

Although most objects in GASPS remain unresolved by 
PACS, five targets (identified in Appendix A) in Taurus were 
found to have clearly extended [OI]63 pm emission along 
known optical jets (Podio et al. 2012). Two of these also have 
broad line profiles in the centre. Figure 7 compares the spectrum 
of one example (RW Aur) with the unresolved line from AATau 
(in red), revealing a prominent red-shifted wing in RW Aur ex- 
tending as much as +200 kms -1 from the stellar velocity. By 
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microns 

Fig. 7. — Spectrum of [OI]63 fim and H 2 0 from the jet source RW Aur ( blue 
histogram), compared with AATau ( red histogram, scaled up by a factor of 3 for 
easier comparison). These spectra are continuum-subtracted, and only the cen- 
tral spaxel is shown. AATau is unresolved whereas RW Aur has a prominent red- 
shifted wing in [01]. The velocities of the three components which dominate the 
[OI]6300 A line are shown by the green tick marks, at the stellar velocity (he- 
liocentric velocity +23 km s -1 ) and at —190 and +100 km s -1 (Hartigan et al. 
1995). Also shown is the wavelength of the H 2 0 line, at the stellar velocity. 


contrast, AA Tau has emission centred at the stellar velocity, 
with a fitted linewidth of 93 kms -1 (FWHM) — similar to the 
measuredPACS resolution of 88 kms _1 at this wavelength 
(see PACS User Manual). The optical [OIJ6300 A line from 
RW Aur is known to originate from highly-excited gas in a 
jet of length a few arcsec (Melnikov et al. 2009), and the line- 
profile is dominated by three components (marked in Fig. 7), 
two at high velocities (+100 and —190 kms -1 ), and one at the 
stellar velocity (Hartigan et al. 1995). The brighter red-shifted 
optical component corresponds with the [OI]63 //m wing, sug- 
gesting this is also from the shocked outflow gas (see § 4.2). 
However, the FIR line profile is dominated by emission centred 
approximately on the star, whereas this velocity component in 
the optical line is relatively weak (Hartigan et al. 1995). This 
low-velocity gas may be from the disk or disk wind (see above). 

In the GASPS data we have also identified five other objects 
with evidence of either broadened lines or spatially-extended 
[OI]63 //m emission: HE Tau and XZ Tau (in the same PACS 
field), DO Tau, UZ Tau, and DK Cha (for the latter source, see 
van Kempen et al. [2010]). All targets resolved in [OI] are iden- 
tified in Appendix A by the note ‘ext.OI’, and it is likely that 
[OI] emission is dominated by outflow gas in these cases. 

5.3.2. Objects with Uncertain Origin of [OIJ63 fim 

In addition to the 10 resolved objects above, a further ~17 
objects (noted as ‘jet’ sources in Appendix A) were identified as 
having published evidence of a high-velocity jet or outflow 
(Kenyon et al. 2008; Podio et al. 2012; Howard et al. 2013, 
in preparation). These are sources with a jet imaged in optical 
lines, a high velocity molecular outflow, or a broad (>50 km s -1 ), 
typically blue-shifted, emission line profile in [O I] 6300 A (see, 
e.g., Hartigan et al. [1995]). Three of these were HAeBe stars 
(HD 163296, MWC480, and HD 100546), leaving 14 T Tauri 
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’jet’ sources. As noted in § 4.2, the survey sensitivity should 
allow us to detect [OI]63 jim emission from outflows shocks 
with mass loss rates M > 3 x 10 -9 M 0 /yr. Estimates from 
Hartigan et al. (1995) suggest that 4 of the jet sources in 
the Taurus sample have mass loss rates exceeding this limit 
(indicated in Appendix A by the note ‘high M’). However, 
their [OI]63 /im emission is neither spatially nor spectrally 
extended in the PACS data. This suggests that the outflow 
shock contribution may be small compared with the low- 
velocity gas; and, the origin of the unresolved [OI]63 /im 
emission in these remaining ‘jet’ sources is not clear from 
the pACS data alone. 

5.4. [CII] Emission 

The detection rate of [CII] 157 /im in the survey was relative- 
ly low. For example, neither of the disks around AA Tau and 
HD 135344 were seen, yet both of these are among the most 
massive disks in the survey (total masses of ~10 -2 M 0 ), with 
relatively rich spectra at other wavelengths. Woitke et al. (2010) 
predicted that the [CII] disk detection rate for Herschel/GASPS, 
assuming a wide range of grid parameters, should be 10-55%, 
and would be highly dependent on the UV excess (§ 4.1.2). 
Table 5 indicates a detection rate at the low end of this range: 
the brightest 44% of [OI]-detected objects were targetted for 
[CII] and of those, only ~26% were detected in [CII]. This 
might indicate that the low-UV models are more applicable 
to the sample. However, this is not supported by the [OI] detec- 
tion rates, which are more consistent with moderate UV ex- 
cesses (§ 5.5.1). Further investigation of this discrepancy is 
warranted. 

A few objects showed extended [CII], or evidence of emis- 
sion from the chop reference position, but in general problems 
from such confusion were limited (§ 2.1.4). There were, how- 
ever, clear cases of both high- and low-mass objects with [CII] 



microns 

Fig. 8. — Spectra of [CII] 157 from UY Aur in the central spaxel ( blue histo- 
gram), and average of the 8 adjacent and 16 outer spaxels ( red and green his- 
tograms, both scaled up by a factor of 6 for clarity). The emission in both line 
and continuum is centrally -peaked and consistent with the PSF response and the 
line-to-continuum ratio is similar in the central and first ring of spaxels, indi- 
cating that neither the line nor continuum are spatially extended compared with 
the beam. 
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emission centred on the star, examples being UY Aur and 
HD 100546. Figure 8 compares the spectra of UY Aur from 
the central spaxel with the middle ring of 8 and outer 16 spaxels 
in the PACS IFU. Both line and continuum are centrally peaked, 
with average fluxes consistent with the instrumental PSF size 
(11") at 157 /im. However, published coronograph images 
shows that these objects also have scattered light extending over 
5-10" (Hioki et al. 2007; Ardila et al. 2007) with a complex 
scattering morphology. This is larger than typical disk sizes 
and suggests emission may be from a compact envelope. The 
origin of [CII] line in these objects and whether it arises from 
the disk, compact envelope, or unresolved outflow is under 
investigation. But it suggests that the [CII] detection rate from 
the disks themselves might be even lower than indicated in 
Table 5. 

5.5. Effect of System Parameters on Line Detectability 

The GASPS target list (Appendix A) comprises a rich sample 
of Class II — III objects in the ~0.3-30 Myr age range, and the 
survey detects [OI]63 //m from half of the targets observed. In 
the following, we investigate preliminary trends in line detect- 
ability vs. other directly-observed parameters. Results from 
GASPS papers on the individual associations (both published 
and in preparation) are combined to look at overall detection 
statistics. A minimum detection limit is 3cr and, although these 
were not all reduced with the same version of HIPE, the criteria 
for detection/non-detection is considered robust in this study. A 
more detailed investigation of correlations of line fluxes using 
systematically-calibrated data obtained from the same software 
version is left for a later paper. 

5.5.7. Disk Dust Mass and [OI] Detections 

The probability of [OI] detection in the GASPS survey is a 
strong function of the disk dust mass, M d . This is illustrated in 
Figure 9, as a histogram of the detection rates as a function of 
distance-normalised dust mass, M' d = M rf .(140/79) 2 , where D 
is the distance in pc. In this figure, we normalised the mass to 
an equivalent object giving the same flux at the distance of 
Taurus. We almost always detect the [OI]63 /im line when 
M' d reaches a threshold of >10 -5 M 0 : 84 ± 10% of targets 
were detected above this mass (where the uncertainty is the sta- 
tistical error). This is comparable with the mass detection limit 
of sub-mm continuum surveys (e.g., Andrews & Williams 
[2005]), assuming a standard mass opacity, k v . Assuming also 
an ISM-like gas/dust ratio can be used for all disks, this implies 
a total mass detection threshold for [OI]63 /tm of ~10 _li M 0 . If 
the gas-to-dust ratio is more typically 10 x lower (as has been 
suggested for TW Hya), then the [OI] observations are detecting 
disks with total masses >10 -1 M 0 . As noted in § 5.3.1, some 
[OI]63 /tm emission can be from outflows; on the plot we indi- 
cate in yellow the targets with spatially or spectrally-resolved 
[OI]. Additional targets ( shown in green shading) are those with 



Fig. 9. — Distribution of normalised dust mass in the GASPS sample, illus- 
trating the mass threshold for [OI]63 fim detection of a disk at the fiducial dis- 
tance of 140 pc. Targets with [OI] detections are shown in shaded colours. The 
dust masses (in Solar units ) are mostly based on published mm continuum ob- 
servations, normalised to the distance of Taurus (140 pc), and assume a standard 
dust mass opacity (see text). Yellow shading indicates objects with extended 
energetic outflows, where the [OI] line is spatially or spectrally extended 
and includes some contribution from the jet (see § 5.3.1 and 5.3.2). Green shad- 
ing indicates objects with evidence of an optical jet, but without spatially or 
spectrally-resolved [OI]63 fim emission; the contribution to the line from the 
jet in these objects is unclear. 


published evidence of a high-velocity jet, although the contri- 
bution of this to the [OIJ63 fim emission flux is unclear 
(see § 5.3.2). 

The 84% detection rate for systems of M' d > 10 -5 M Q drops 
to 32 ± 12% for 1(T 6 < M' d < 10~ 5 M 0 . Woitke et al. (2010) 
constructed a large grid of disk models covering a parameter 
space similar to the that of the GASPS sample, and predicted 
overall [OI]63 //in detection rates of 51-70% for disks with dust 
masses of 10' ' — 10 — 3 M e in systems with a high UV excess, 
and 17-30% for this mass range in the case of low UV. Restrict- 
ing the model grid of Woitke et al. (2010) to include only disks 
of dust mass >10 -5 M G , we find that such [OI] detection rates 
are achieved for moderate UV excesses (0.01 < / uv < 

0.1) and disk flaring (1.0 < (3 < 1.2). This suggests these 
ranges are typical of most systems in the survey. 

There are some notable exceptions to the mass detection 
threshold, where we detected [OI] in systems with M' d < 3x 
10“ 6 M 0 : 

1. HD 172555, an unusual warm debris system with no evi- 
dence of molecular gas in mm lines, but with some indication 
that [OI] may be secondary gas released in collisions (Riviere- 
Marichalar et al. 2012b). 

2. ET Cha, an apparently compact disk in the relatively old r/ 
Cha association (Woitke et al. 2011). 
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Teff 

Fig. 10. — Disk dust mass of the GASPS sample, normalised to a distance of 
140 pc ( units ofM a ), plotted as a function of the stellar effective temperature (in 
K). This illustrates the line detectability of a disk of a given dust mass at the 
distance of the Taurus star forming region. [OI]63 //m detections are shown as 
filled circles, and open circles depict [01] upper limits. An additional red dot 
indicates systems which were detected in [CII], 


3. J 130521 . 6-773810, although the classification of this tar- 
get in Chall is uncertain. 

4.51 Oph, a warm compact disk with notable hot and com- 
pact molecular gas component (Thi et al., submitted). 

5. HD141569, a diffuse disk with spiral structure around a 
HAeBe star (Clampin et al. 2003). 

The number of disks with published dust masses as low as 
10 -8 M 0 is relatively small (only ~10 in GASPS have 
measured values), and further mm-wavelength measurements 
of such disks would be interesting to improve the statistics. At 
the opposite extreme, three relatively massive disks ( M' d > 
10 -4 M 0 ) have no evidence of [OI]63 //in: GO Tau, V836 
Tau, and TWA03. Woitke et al. (2010) indicate that disks of this 


mass which have low flaring (/3 < 1.0) can have [OI]63 /im 
fluxes too low to be detected by GASPS. 

5.5.2. Dependence on Spectral Type 

It is already clear from Table 5 that [OI]63 /im is significantly 
easier to detect around HAeBe stars than T Tauri stars. Is this 
simply because HAeBe disks in the sample are more massive 
and the detection threshold is more commonly reached? 
Figure 10 shows the distribution of normalised disk dust masses 
{M' d ) in GASPS as a function of stellar T eff . Systems detected in 
[OI] and [CII] are indicated by the filled black and red symbols 
respectively. This shows that both early and late-type stars have 
a similar range of disk dust masses in this sample. As noted 
above, the [OI] detection rate is high for disks with M' d > 
10 -5 Mq, and Figure 10 shows that this is independent of 
T eff for T eff > 4000 K. However, approximately half of the 
low-luminosity stars (T eff < 4000 K, or M type) with M' d in 
the range 10 -5 to 10 -4 M Q were not detected. Clearly the spec- 
tral type has some effect on the [OI] line emission threshold for 
the lowest-luminosity stars. 

In the case of the [CII] 157 line, the detectability in Figure 10 
seems to be independent of the spectral type and disk mass, with 
[CII] detections (filled red symbols) broadly distributed over the 
M' d -T ei{ parameter space. Unlike [OI], there is no clear thresh- 
old with disk mass, or an increase in detection rate among 
HAeBe stars. If most [CII] arises from a compact envelope rath- 
er than the disk, this suggests that such gas may be retained 
around these stars independent of the mass of the inner disk 
or stellar type. 

5.5.3. Other Observational Parameters: Binarity, Ha., and 
X-ray Luminosity 

The histograms in Figures 1 1 and 12 show the detection sta- 
tistics for targets searched in [OI]63 pm with published X-ray 



Fig. 11. — Flistogram of the distribution of (a) X-ray luminosity (erg s 1 ) and (b) Fla EW (A) of the GASPS sample observed in [OI]63 /im. The [OI] detections are 
shaded in blue, with the yellow shading indicating those with an extended [OI] component. There is no clear dependence of [OI] detectability on X-ray luminosity, but 
detection rates are higher for larger Ha EW. 
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Separation (AU) 


Fig. 12. — Histogram of the distribution of binary separation (in AU) in the 
sample observed in [OI]63 fim, with line detections shaded in blue. Binary stars 
with known extended jet contribution to the [OI]63 /im flux are shaded further in 
yellow. 

luminosity, Ha equivalent width (EW) and binary separation. 
The numbers of targets with [01] detections are shaded. Those 
which have additional extended [OI]63 /im emission from a jet 
are shaded yellow. Figure 11a shows that line emission is de- 
tected in systems covering the full range of X-ray luminosity in 
the survey, with no clear trend of increased detectability for 
higher X-ray fluxes. The Ha EW used in Figure 11 b is linked 
with the accretion rate, although later-type K stars may have 
significant chromospheric contribution and the accretion lumi- 
nosity may be lower than Figure 11 b might suggest. But there is 
a trend of increasing [01] detection probability for higher accre- 
tion rates, including all stars observed. The detection fraction is 
70% for EW >30 A (or 67% excluding the stars with extended 
jet emission) compared with only 29% for those with lower EW. 
A systematic derivation of accretion luminosity and line flux 


over the whole survey would be interesting to study further cor- 
relations between the [OI]63 /im flux and accretion rates. 

Figure 12 illustrates the [01] detection rates distributed over 
binary separation. For hierarchical multiples we have used the 
separation of the widest component within the PACS beam. 
There is marginal evidence for a drop in detection rates in multi- 
ples of separation <300 AU, from 64% for the wider binaries to 
40% for the closer systems (with statistical errors of — 10%). By 
comparison, the [OI] detection fraction of single stars in the 
sample was 47%. This would suggest that most [OI] emission 
arises from radii of <300 AU — similar to the [OI] -emitting 
region suggested by the models in Figure 1. By comparison, 
samples of T Tauri stars observed in mm dust indicate that 
M d typically drops by a factor of 5 for binary separations of 
<300 AU (Harris et al. 2012). 

5.5.4. Detection Rates in Different Associations and 
Dependence on Age 

In Table 6 we give the [OI]63 /tm detection rates for the dif- 
ferent associations observed in GASPS, for those with and with- 
out published dust masses (for masses based on continuum 
photometry, although only two disks were detected which had 
no published estimates). As noted above, the required dust mass 
( M ' d ) for [OI]63 /tm to be detected at a distance of 140 pc is 
~10 -5 M q , or a total (gas+dust) mass of 1 Mj upiter , assuming 
an ISM gas-to-dust ratio. The detection rates are somewhat de- 
pendent on spectral type (M stars have a higher disk mass 
threshold for detection — see Fig. 10), binarity (lower for separa- 
tions <300 AU), and accretion rates (more for H a > 30 A). 
However, if we use [OI]63 /im detections as a proxy for 
1 M Jupiter disks at the fiducial 140 pc distance, it is possible 
to compare the detection rates in different associations, modulo 
the distances, and assume we are sampling most of the brighter 
detectable disks. Each of the three intermediate-age (5-10 Myr) 
associations (77 Cha, TWA and upper Sco) has ~2-3 such disks. 
TWA is significantly closer, and only one of these would be 


TABLE 6 

Detection Statistics of [OI]63 /im for the Associations in GASPS, With and Without Known Disk 

Masses 

M d No M / Total 


Association Detected Observed Detected Observed Detected Observed 


Taurus 44 56 1 17 45 73 

Cha II 7 17 1 2 8 19 

r] Cha 2 7 0 6 2 13 

TWHya 3 5 0 3 3 8 

Upper Sco 2 7 0 12 2 19 

.1 Pic 1 4 0 1 2 b 6 

Tuc Hor 0 2 0 0 0 2 

HAeBe stars 20 25“ 0 0 20 20 


“No disk masses were available from the literature. 
b Includes detection of /3 Pic itself (Brandeker et al. [2011]). 

“Five nondetected systems are those classed as possible debris disks. 
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detected by GASPS if moved to the fiducial distance of 140 pc. 
The total number of stars in the upper Sco group is ~200 
(Carpenter et al. 2009), with 20-50 in the two other associations 
in this age range (Torres et al. 2008), giving a gas-detected disk 
fraction of 1-7% at 5-10 Myr. For the older systems, there are 
two unusual [OI]-emitting disks in the 10-20 Myr f3 Pic moving 
group but neither would have been detected at 140 pc, and no 
disks were detected in [OI] in the 30-Myr old Tuc Hor associa- 
tion. Each of these contains ~40 stars (Torres et al. 2008). For 
the two younger groups, the detection rate in GASPS was 
~50%. However, the total number of stars is estimated to be 
~250 in Taunts and ~48 in Cha II (Rebull et al. 2010; Alcala 
et al. 2008), giving massive disk fractions of ~18 and 17%, re- 
spectively. For the HAeBe stars in the sample (excluding debris 
disks), the derived ages are mostly in the range 2-10 Myr 
(Montesinos et al. 2009; Meeus et al. 2012), and the [OI] frac- 
tion is ~100%, much higher than the equivalent-aged FGKM 
star. While not a statistically-complete sample of AeBe stars 
in this age range, it suggests either that these more luminous 
stars are more commonly able to retain disks, or that their ages 
are overestimated. Overall, for T Tauri stars, the fraction with 
Jupiter-mass, gas-rich disks is ~18% at ages of 0.3-4 Myr, 
1-7% at 5-10 Myr, and none are detected beyond 10 Myr 
age. It is unclear why some particular stars can retain these 
gas-rich disks for up to 10 Myr, and whether planets formed 
in such long-lived disks would be different from those around 
other stars. 

6. CONCLUSIONS 

In this article we describe the GASPS far-infrared survey of 
gas and dust in young stellar systems. This Herschel Key project 
observed selected lines and continuum with the PACS instru- 
ment, targetting 250 young SED-Class II— III systems, Herbig 
AeBe stars, and young debris disks. The far-IR line emission 
can arise from the hot surface of gas disks around 30-300 AU, 
high-velocity jet shocks, disk winds and/or compact remnant 
envelopes. We outline the models used to predict disk line 
strengths and their dependence on disk parameters, and look 
at the possible emission from these other mechanisms, in par- 
ticular from extended outflow jets in a few objects. The master 
list of targets with basic system parameters is presented, and we 
also indicate in this list which objects have detections in the 
primary GASPS lines. 

The FIR lines are identified and we show preliminary line 
detection statistics, referring to the published papers which an- 
alyse individual sources and associations in more detail. Rich 


spectra were seen in a number of sources, including fine- 
structure atomic lines of [OI] at 63 and 1 45//m, and [CII] at 
157 ^m, as well as molecules including high-J transitions of 
CO, H 2 0, OH and CH + . In most systems, [OI]63 //m was 
the brightest line by a factor of ~10, and is associated in most 
cases with disk emission. Line and continuum emission was, in 
all but 10 systems, both spatially and spectrally unresolved and 
centred on the star. For the extended sources, emission is 
thought to arise from jet shocks, along with a disk contribution. 
The [CII] 157 line was significantly weaker than [OI], resulting 
in a relatively low detection rate. However, in a small number of 
objects unresolved emission was found centred on the star, and 
may arise from a compact remnant envelope. 

49% of all targets observed were detected in [OI]63 /im, with 
an 84% detection rate for those having disk dust masses of 
>10 -5 M 0 . Comparison with statistics from a grid of model 
implies that most systems have moderate UV excesses and disk 
flaring. Assuming an ISM gas-to-dust ratio of 100 and typical 
mm-wavelength mass opacity, this represents an [OI] detection 
threshold for the total disk mass of ~1 M Jupiter . Going against 
this trend, we find five unusually low disk mass systems with 
[OI]63 /im emission, and a few high-mass systems which re- 
mained undetected. The [OI]63 //m detection rates were inde- 
pendent of X-ray luminosity, but there was evidence of a 
decreased rate in binaries of separation <300 AU, in stars with 
Ha EW <30 A, and for M-type stars. Based on [OI]63 //m line 
detections, the results show that ~18% of stars in each associa- 
tion can retain gas-rich disks of minimum mass lO"' 5 M 0 for up 
to 4 Myr, a few % of stars keep these disks for ~5-10 Myr, but 
none are detected beyond 10-20 Myr. 
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2009-256513), CNES, and PNPS of CNRS/INSU, France for 
support.PW, IK, and WFT acknowledge funding from the EU 
FP7-2011 under Grant Agreement nr. 284405. WFT acknowl- 
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LAC (Spain). This development has been supported by the fund- 
ing agencies BMVIT (Austria), ESA-PRODEX (Belgium), 
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